Interplay of friction and noise and enhancement of disoriented chiral condensate 
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Using the Langevin equation for the hnear a model, we have investigated the effect of friction and 
noise on the possible disoriented chiral condensate formation. Friction and noise are supposed to 
suppress longwavelength oscillations and growth of disoriented chiral condensate domains. Details 
simulation shows that for heavy ion collisions, interplay of friction and noise occur in such a manner 
that formation of disoriented chiral condensate domains are enhanced. 
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In recent years there is much excitement about the possibility of formation of disoriented chiral condensate (DCC) . 
The idea was suggested by Rajagopal and Wilczek y] and also by Bjorken and others In hadron-hadron or in 
[ heavy ion collisions, a macroscopic region of space-time may be created within which the chiral order parameter is not 
CsJ . oriented in the same direction in the internal 0(4) = SU{2) x SU{2) space as the ordinary vacuum. The misaligned 
^ " condensate has the same quark content and quantum numbers as do pions and essentially constitute a classical pion 
Q , field. The system will finally relaxes to the true vacuum and in the process can emit coherent pions. Possibility 
of producing classical pion fields in heavy ion collisions had been discussed earlier by Anslem DCC formation 
in hadronic or in heavy ion collisions can lead to the spectacular events that some portion of the detector will be 
dominated by charged pions or by neutral pions only. In contrast, in a general event all the three pions {tt~^,tt~ and 
7r°) will be equally well produced. This may then be the natural explanation of the so called Centauro events 

Microscopic physics governing DCC phenomena is not well known. It is in the regime of non-perturbative QCD as 
(y-^ ' well as nonlinear phenomena, theoretical understanding of both of which are limited. One thus uses some effective 
field theory like linear a model with various approximations to simulate the chiral phase transition ||^-|^ . In the linear 
sigma model chiral degrees of freedom are described by the the real 0(4) field <I> = (tr, tt). Because of the isomorphism 
between the groups 0(4) and SU{2) x SU{2), the later being the appropriate group for two flavour QCD, linear sigma 
model can effectively model the low energy dynamic of QCD Explicit simulation with linear sigma model, indicate 
that DCC depends critically on the initial field configurations. With quench like initial condition DCC domains of 
4-5 fm in size can form Q. Initial conditions other than quench lead to much smaller domain size. Quench scenario 
assume that the effective potential governing the evolution of long wave length modes immediately after the phase 
transition at Tc turns to classical one at zero temperature. It can happen only in case of very rapid cooling and 
O ^ expansion of the fireball. In heavy ion collisions quench like initial conditions are unlikely. 
D , Very recently, effect of external media on possible DCC is being investigated |p|-p^. Indeed, in heavy ion collisions, 
i-C ' even if some region is created where chiral symmetry is restored, that region will be continually interacting with 
^ , surrounding medium (mostly pions). The surrounding medium can be conveniently modelled by a white noise source 
and one can use Langevin equation for linear a model to simulate the DCC formation under the infiuence of external 
media. Recently it has been shown that in 0* model, hard modes can be integrated out to obtain a Langevin type of 
equations for the soft modes Biro and Greiner ||l^ using a Langevin equation for the linear a model, investigated 
the interplay of friction and white noise on the evolution and stability of zero mode pion fields. In general friction 
and noise reduces the amplification of zero modes. But in some trajectories, large amplification may occur | p^ . We 
also obtain similar results fl^ . 

It is popular expectation that friction and noise will reduce the DCC formation probability. The expectation was 
found to be true in one dimensional calculations with zero modes only ]To| , p^ . To see how far this expectation is 
valid when higher modes are included, in the present letter we have investigated DCC formation in 3-f 1 dimension. 
Two scenarios were considered. In scenario I, we use the equation of motion of linear a model with quenched initial 
condition to simulate DCC formation without friction and noise. In scenario II, we solve the Langevin equation for 
linear a model, with the same quenched initial condition to simulate DCC formation in presence of friction and noise. 
Comparative study of these two scenarios reveal that certainly the noise and friction affect DCC formation probability. 
However contrary to one dimensional calculations, large DCC formation is seen in scenario II rather than in scenario 
I. Indeed, it will be shown that large DCC domains are formed in scenario II only, not in scenario I. This effect is 
particular to heavy ion collisions as explained below. 

The Langevin equation for linear a model, at temperature T can be written as, 

1^ + (7 + + A(*= - f, - T'/im - ar...y.Y) (1) 

where = (tr, tti, 7r2, tts). t is the proper time and Y is the rapidity, the appropriate coordinates to describe heavy 
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ion collisions, t^o-.tt and C,a,TT are the friction coefficients and the white noise for the a and tt fields. We note that if the 
friction and the noise terms are dropped from eq.|l| the resulting equation is for the scenario I, i.e. just the equation 
of motion of linear cr model. 

The noise source C, and the friction -q are not independent. They are related by the fluctuations-dissipation relation. 
We use white noise source with zero average and correlation as demanded from the fluctuation-dissipation relation, 

<Ca{T,x,y,Y)>^Q (2a) 
< Ca(T, X, V, Y)Ur'.x\ v\ Y') >= 2Tf^-5{T - t')5{x - x')5{y - v')5{Y - Y')8ab (2b) 



where a, h corresponds to tt or a fields. A few words are necessary about the use of temperature in the noise term. 
We are approximating DCC, which is a non-equilibrium phenomena as a equilibrium one. Such an approximation is 
valid when the system is not far from equilibrium. Indeed, fluctuation-dissipation theorem is valid for such a system 
only. 

Friction coefficients is an important ingredient for the Langevin equation. In our earlier study, we have assumed that 
both a and tt evolve under the influence of a common friction [|l5|,^ . In the present work we treat them separately. 
a and tt fields evolve under the influence of friction, appropriate for them. Friction coefficients for a and tt have been 
calculated by Rischke p^, 

^ ....^ l-exp{-m^/T) 1 

^ N ' A-KmlV ml 1 - exp{-ml/2mnT) exp{ml - 2ml)/2m^T) - 1 ^ ' 




A\f^ N-l 4m2 
''-^^^^^^^ ■ \^ -^^^^-^ (3^) 



Solution of eq.|i| require initial fields configurations. They were distributed according to a random Gaussian with, 

<a>={l- ,f{r))U (4a) 

< TTj >= (4b) 

<a^>-<a >2=< Trf > - < tt, /J/6/(r) (4c) 

< a < TTi (4d) 

< >=< 4/2/6/(r) (4e) 



The interpolation function 

fir) = [1 + exp(r - ro)/T)]-' (5) 

separates the central region from the rest of the system. We have taken ro=6.4 fm and r=.7 fm. The initial field 
configurations corresponds to quench like condition ^ and as told in the beginning arc unlikely to be obtained in a 
heavy ion collisions. We still use it as they are the most favourable initial conditions to produce DCC like phenomena. 
In the simulation results presented below, the initial time and temperature were assumed to be Ti — lim/c and Ti=200 
MeV. Effect of expansion was included through the cooling law, 

T{t)=tA (6) 

T 

which is rather fast for heavy ion collisions. However, we choose to use it as DCC formation probability increases if 
the system cools rapidly. 

With the initial conditions and the cooling law most appropriate for DCC formation, we have solved the Langevin 
equation (scenario II) on a 32'^ lattice, with lattice spacing of a=l fm, using a time step of a/10 fm/c. We also 
use periodic boundary conditions. The equation of motion for linear a model (scenario I) was solved similarly, with 
identical initial conditions. We have continued the evolution of the fields for about 10 fm/c. 

We define a correlation function at rapidity F as 
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where the sum is taken over those grid points i and j such that the distance between i and j is r. In fig.l, we have 
compared the correlation function in scenario I and II. Initially there is no correlation length beyond the lattice 
spacing of 1 fm. Correlations starts to develop at later time. It increases for about 7fm/c, then decreases again. 
Interestingly, larger correlation length is obtained in the scenario II, than in scenario I. Thus at 7 fm/c, correlation 
length in scenario I is only 2 fm, while that in scenario II is 4 fm. Increased correlation in scenario II is contrary 
to popular expectation that friction and noise will reduce correlation. 

Enhancement of correlation length in scenario II, with friction and noise is corroborated in the pion field distribution 
also. In fig. 2, we have shown the xy contour plot of the tt2 component, at rapidity Y=0. Field distribution at Ti=l 
fm/c and after 7 fm/c of evolution are shown. Initially there is no correlation. Domain like structure is seen both in 
scenario I and II, after 7fm/c evolution . The positive and negative components of the tt2 separates out. Here again, 
much larger domains are formed in scenario II than in scenario I. It may be noted that domain like structure seen 
in one of the component of tt field do not necessarily convert into physical domains. They are indication of larger 
correlation length only. Physical domain should contain either charged or neutral pion only. Thus in physical domain 
neutral to total pion ratio should differ considerably from the isospin symmetric value of 1/3. 

Assuming that the pion density is proportional to the amplitude's square, in fig. 3, we have shown the contour plot 
of the neutral to total pion ratio, at rapidity Y = 0. 

i?3(x,y,F,r) = -^^|^ (8) 

Very small or large value of the ratio over an extended spatial zone will be definite indication of domain formation. 
As expected, initially there is no domain like structure. In scenario I, we donot find any large domain like structure 
in the ratio R3 even after 7 fm/c evolution. Thus while tt2 component of the tt field shows domain like structure at 
7 fm/c, in terms of physical pious there is no domain structure in scenario I. Large DCC domains are not formed in 
linear a model, even with quenched initial condition and fast cooling law. In scenario II however we can find some 
extended zone with large or very small value of the ratio i?3 . Thus physical pious evolve into domain like structure 
in scenario II, with noise and friction, rather than in scenario I. 

Present simulation indicate that friction and noise do indeed enhance DCC domain formation possibility. The 
simulation results are contrary to expectation from one dimensional calculations. As such noise and friction are 
supposed to reduce the amplification of long wavelength modes. We believe that present results are particular to 
heavy ion collisions where proper time and rapidity are the most appropriate coordinate systems. In this coordinate 
system, correlation of the noise decreases with (proper) time (1/r dependence, see eqj|). Physically the source of 
noise i.e. the medium surrounding the zone where chiral symmetry is restored, fly away with time, reducing the 
correlation. However, the friction coefficients 770-, tt remains more or less same in the temperature range considered. 
Thus at later time evolution of the fields are determined mainly by the friction. We also not that ryo- is considerably 
large (~ 3/m~^). Then once the trajectory enters into unstable regime, large friction opposes its tendency to come 
out of the instability. Friction forces the trajectory to remain in the unstable regime for longer duration. Indeed, in 
one dimension, we have obtained similar result 3d simulation confirms our results in one dimension. 

Though noise and friction enhances the DCC formation probability, it may not be easy to detect it. In fig. 4, we 
have shown the rapidity distribution of the neutral to pion ratio, 

_ jTrldxdy 



at different time intervals. Through out the range of rapidity, the ratio fluctuate about its isospin symmetric value 
of 1/3. As expected fluctuations are larger in scenario II than in scenario I. However, if we remember that the pious 
that will be detected are integrated over time, then it is obvious that the fluctuations will be considerably less. Thus 
even in event-by-event analysis, it will be difficult to tell about DCC formation from rapidity distribution of pious 
only. Much more study is needed to resolve the issue. 

In summary, we have considered disoriented chiral condensate domain formation in two scenarios, one without 
noise and friction (scenario I) and the other with noise and friction (scenario II). In scenario I, equation of motion 
for the linear a model fields and in scenario II, the Langevin equations for linear a model were solved. Using the 
most ideal conditions for DCC formation, i.e. quenched initial condition and fast cooling, we have evolved the fields 
for 10 fm/c. It was seen that in both the scenario, tttt correlation increases with time till 7 fm/c, it then decreases. 
Larger correlation is obtained in scenario II than in scenario I. Evidence of increased correlation is also obtained from 
the contour plot of pion field. 7r2 component of the pion field shows domain like structure. Positive and negative 
component of fields separates out at late time. Here again, domain like structure is more prominent in scenario II 
than in scenario I. Contour plot of the neutral to total pion ratio suggest that though domain like structure is seen 
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in one component of the tt field, there is no (large) physical domain formation in scenario I. On the contrary, large 
physical domains are seen to be formed in scenario II after 7 fm/c of evolution. We have also studied the rapidity 
distribution of the neutral to pion ratio in both the scenarios. Throughout the rapidity range, the ratio fluctuate 
about the isospin symmetric value of 1/3, fluctuations being more in scenario II than in scenario I. However, it was 
also noted that fluctuations will be considerably less when integrated over time. Rapidity distribution of neutral to 
total pion ratio may not be able to single out the DCC events. 
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FIG. 1. Evolution of the pion correlation fucntion with time, in scenario I and II. 
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(a) 7t at X = 1 fm 
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(b) 71 at X =7 fm, scenari I 
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(c) 71 at X =7 fm, scenari II 
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FIG. 2. Contour plot of the tt2 component of the pion field at rapidity Y=0, at initial time Ti — l fm/c and after evolution 
for 7 fm/c. 
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(a) R3 at X = 1 fm 
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FIG. 3. Contour plot of the neutral to total pion ratio, at rapidity Y=0, at initial time ri=l fm/c and after evolution for 7 
fm/c. 
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FIG. 4. Rapidity distribution of neutral to total pion ratio, in scenario I and II. The ratio at different time intervals are 
shown. 
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